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ABSTFACT 

Accurate  prediction  of  variations  in  the  height  of 
atmospheric  inversions  are  required  for  optimum  utilization 
of  modern  weaponry.  A  marine  atmospheric  boundary  layer 
model  was  used  to  explore  the  influence  of  land  or 
predictions  of  the  inversion  heights.  Data  for  this  evalua- 
tion were  cfatained  from  rawinsonde  observations  taken  at 
Vandenfcerg  AFB,  Califcrnia  for  each  season.  Vertical  veloc- 
ities fcr  these  data  were  obtained  by  vertical  integration 
cf  the  mcisture  budget  equation  (Q-method) .  Results  of  the 
model  output  were  compared  graphically  and  by  use  of  RMS 
error  statistics.  These  results  were  very  poor  and  revealed 
that  strong  land  influences  on  diurnal  variations  yield 
inversion  height  fluctuations  which  the  model  cannct  repro- 
duce, further  development  of  the  model  is  necessary  to 
overcome    problems  caused  by   the   land   influence. 
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Historically ,  advanced  detection  and  surprise  attacks 
have  teen  an  area  cf  interest  to  tha  military  community. 
Knowledge  about  the  effects  of  atmospheric  boundary  layer 
(A3L)  inversions  on  electromagnetic  signals  is  well  known 
and  documented.  Advance  knowledge  of  inversion  changes  can 
have  a  direct  effect  on  military  planning  conducted  by  ail 
services. 

Kith  technological  advances,  our  weaponry  has  become 
highly  sophisticated  (e.g.,  the  development  of  "smart" 
bombs)  and  dependent  on  electromagnetic  propagation  for 
guidance  ccntrol.  Tfce  guidance  components  for  these  weapons 
deperd  en  many  different  wavelengths,  ranging  from  the 
visible  to  the  microwave,  dependent  on  the  specific  use  of 
the   weapon.  The    weapons      defend   on      the    transmission      and 

reception  of  these  signals  to  locate  and  follow  targets. 
Difficulties  arise  wten  atmospheric  conditions  are  such  that 
signal  performance  is  distorted  by  refraction.  Only  these 
areas  which  affect  propagation  of  radar  signals  will  be 
discussed  in  this  thesis.  Other  propagation  types,  such  as 
optical,  and  the  atmospheric  parameters  affecting  them  will 
rot   be    addressed. 

Signal  propagation  is  dependent  on  the  index  of  refrac- 
tion which  is  a  function  of  temperature  and  humidity. 
Eefracticn  of  a  sigral  occurs  as  that  signal  moves  through 
vertical  gradients  of  temperature  and  humidity.  Atmospheric 
conditions  which  enhance  refraction  are  inversions,  which 
cause  strong  gradients  of  temperature  and  humidity,  or  rapid 
changes  cf  humidity  near  the  surface.  Occurrence,  height 
and  intensity  of  inversions  are  strongly  dependent  on  local 
and      synoptic      atmospheric    conditions.         The      importance      of 
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predicting    propagation    is   ofcvicus   but   to   do   sc   wi 
first      reguirss    accurate      prediction      of    the 

r*  p  T*  —  m  rr\  +-  .~  "J-  c  ^  -f"    £  *-  {•"**    1  •n  f       *   "^  Q  T*  £  h  VI  O  U  ^  *    ft]  D  O  "*""*"  "*    T"1  *~ "  ~         *     r  **  V    £ 

prediction  of  changes  in  the  temperature  and  humidity 
profiles,  both  of  which  are  highly  variable  and  can 
completely   change   their   vertical   profile      within    a   matter   of 


hours 


Other    paraseters      such      as      vertical   velocity      and 


surface    v»ind  speeds    Kill   be    discussed    later. 

The    ability    to    forecast    inversions   and    how   they    iray    vary 


directly   influences 


tn< 


ffectiveness      of    modern      weaponry, 


It  is  readily  apparent  that  improvements  in  this  forecasting 
ability  is  cf  paramount  importance  in  making  military  plans 
which  invclve  choosing  the  type  weapon  best  suited  for  a 
particular  mission.  As  long  as  there  are  functions  depen- 
dent en  signal  transmissions  there  will  remain  a  need  to 
increase  the  ability  tc  forecast  the  atmospheric  parameters 
which   contribute    to    propagation. 

The  model  used  tc  investigate  inversion  heights  is  a 
narir.e  atmospheric  boundary  layer  (flABL)  model  which  has 
been  modified  for  use  at  the  Naval  Postgraduate  School 
(NPS)  .  The  factors  incorporated  into  this  model  include 
potential  temperature  and  specific  humidity  within  the 
well-mixed  region  and  their  gradients  at  the  base  of  the 
inversion.  Use  of  these  parameters  make  this  model  appli- 
cable to  studying  inversion  changes.  Actual  discussicn  cf 
the    trcdel    will    be    provided    in    Chapter    II. 

Another  area  where  advance  knowledge  of  ABL  inversion 
changes  t»culd  be  useful  is  in  single-station  forecasting. 
The  effects  of  chances  in  the  inversion  are  found  not  only 
to  effect  electronic  transmissions  but  also  to  indicate 
possible  changes  to  current  weather  conditions.  Therefore, 
any  icprcvements  in  cur  ability  to  predict  these  changes 
will  be  important  tc  land -based  and  sea-based  operations, 
since  improvement  in  this  area  cf  prediction  also  will  be 
reflected    in  the   daily    weather    forecast. 
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This  thesis  completes  an  analysis  on  four  time  periods 
cf  data  obtained  frcm  Vandenberg  AFB,  California.  These 
data  periods  each  cover  a  different  season  and  vary  frcm  two 
to   eight      days    in   length.  In    conducting      this    examination 

vertical  velocities  fcr  Vandenfcerg  AFB  were  calculated  from 
their  ratiinsonde  data.  Each  days  data  were  run  through  the 
NPS  MflEL  prediction  model  to  establish  the  forecast  trend 
for    that   period.  Sesuits    from    the    model's      forecasts   were 

compared  to  data  frci  later  reports  for  correlation  and  to 
determine    accuracy      cf   the    model.  Successive   runs      of    the 

model  were  made  until  each  data  set  has  been  completed  and 
deterninaticn  as  to  how  well  the  model  predicts  inversion 
heights    in    each    24    hcur    forecast. 

One  cf  the  time  periods  of  data  chosen  for  study  is,  by 
choice,  coincident  with  a  data  period  used  by  Gleason 
(1982).  His  study  vias  directed  at  effectiveness  of  the 
model  fcr  mere  of  a  pure  marine  location  and  fcr  this  reason 
data    frci    San      Nicolas   Island    was    chosen.  This    study   uses 

data  from  Vandenterg  AF3  since  it  is  a  coastal  station  and 
will  have  scmewhat  different  surface  conditions  due  to  the 
continental   influence.  See   Fig.         1.1    for     a   geographical 

orientation  of  the  area  of  interest.  Once  the  Vandenberg 
AFB  data  had  been  studied,  the  San  Nicolas  Island  results 
were  used  tc  make  comparisons.  The  results  cf  this  compar- 
ison are  cf  interest  since  they  will  indicate  whether  or  not 
there  are  potential  problem  with  the  tfASL  model  when  used  in 
a  regiie  strongly  influenced  by  land  as  compared  tc  a  regime 
closely   resembling    marine  conditions. 
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figure  1.1    Southern  California  Coastal  Regions. 
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II.    HOD EI    DESCRIPTION 

1.       E1CKGEC0UD 

The  MAEI  model  id  use  at  the  Naval  Postgraduate  School 
(NPS)  is  a  zero-order,  two  layer,  integrated  mixed  layer 
modal  (Davidson  et  al.  1984).  This  is  a  modified  versicn  of 
the  model  developed  by  Stage  and  Businger  (1981).  The  two 
layers  censist  of  th«  wall-mixed,  turbulent  boundary  layer 
and        the        relatively  ncn-turbulent.        free        atmosphere. 

Separating  these  layers  is  the  transition  zone  known  as  the 
inversion.  The  tern,  "zero-crder",  applies  to  this  model 
because  the  inversion  is  considered  a  jump  from  one  value  to 
another  instead  of  having  a  finize  gradient.  Thus,  we 
assume  xhat,  at  the  inversion,  the  equivalent  potential 
temperature  (9  ),  and  the  total  specific  humidity  (q ^  ) 
undergc    zero-thickness    jumps. 

Total  specific  hunidity  and  equivalent  potential  temper- 
ature are  conserved  guantities  in  pseudo-adiabatic  processes 
and,  for  this  reason,  are  assumed  well-mixed  within  the 
boundary   layer.      They  ara   given    by: 


*t    =    %    +    qi 


(2.1) 


where      g      is     the    water      vacor   mixing      ratio   and      q      is 

V  L 

liquid   water   mixing    ratio;     and 


the 


+     L     /  (C     q     ) 
v  p     v 


(2.2) 


where    8    is   the    potential   temperature,      L      is   the    latent    heat 
for    vaporization   of      water,      and   C      is   the      specific    heat   at 


constant    pressure      of   air. 


3   third      fundamental    dependent 
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Figure   2. 1        Typical    potential   temperature   and 
noisture    profiles. 


variable      is     the      in€an   mixed      layer      depth       (h)  .  Typical 

profiles  of  these  parameters  are  shown  in  Fig.  2.1.  As  can 
he  seer  in  Fig.  2.  1,  both  6  and  q.  are  assumed  tc  have 
constant  values  throughout  the  mixed  layer.  Above  this  is 
the  interfacial  regicr,  or  inversion,  where  the  jump  occurs. 
One  cf  the  more  appealing  features  about  this  model  is 
that  it  runs  from  input  parameters  gleaned  frcm  routine 
iteteciclcgical  and  cceancgraphic  data.  These  parameters 
are:  surface      layer      humidity,        temperature      and      winds, 

vertical      gradients    cf      temperature    and      humidity    above      the 
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nixad  layer,  and  sea-surface  temperature  (SST)  .  Readily 
available  rawinscnde  observations  can  be  used  to  determine 
all    input    parameters   except    SST.      SST    values   can    re   obtained 

from    satellite    interpretation    cr   direct   measurements. 

The    predictive      equations    for   the    three      dependent    vari- 
ables   are    given    by: 

,  -1 


dq    /dt     -    h         (      (w«q«)         +     w     Aq     ) 

u  t    u  e       t 


-1 


d9e/dt    =    h        (     [w'8'1.    +    r       -    Rd    +    w    ±Q       )  (2.3) 


e     0 


dh/dt    =    w       +    w 
D  e 


where    d/dt    represents      the    time    rate    of    change      of    the    vari- 
able   following    a   column   of    air.       (w'e')0     and    (w'q^)0    are   the 

surface   turbulent   fluxes   of    6      and    q. ,    while   the    A    terms   are 

e  3t ' 

the    magnitude   of   the    jumps    at    the   inversion.         Other    quanti- 
ties   introduced    here    are:         1) 


w      -    the    entrainraent    rate; 
e 


2) 


-   the   large-scale   vertical   velocity   at   the    inversion 


height  (initialized  using  the  subsidence  rates  described 
later  in  this  chapter);  3)  E_  -  the  rate  of  radiative  heat 
loss    per    unit  area      near    the   cloud    top;      and      4)  Rc    -   the 

rats    of   radiative    heat   gain    per    unit   area    at   the    cloud    base. 
(Haqgerty,     1983) 


E. 


NES    MCDEL    PARTICULARS 


The  use  of  a  well-mixed  layer,  along  with  other  develop- 
ments, irade  modifications  to  the  model  described  by  Stage 
and    Eusinger      (1981)      necessary.  The    NPS      integrated    MABL 

differs  from  others  in  the  specifics  of  surface  flux,  radia- 
tive flux  ard  the  zeic  gradient  assumption.  The  following 
discussions  are  brief  since  the  references  (Davidson  et  al. 
1984)    are    quite    complete. 
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1  •      Atmospheric    structure 

The    profile    cf   a  well-mixed   atmospheric   variable,   $ , 

'  e  tZQ  £31!  ted     a  s 


4>(z)    =4>  +    r      z  z<h 

ms  <pm 

*  (Z)     =    *m*  +    rA  (z~h)  Z>h 

mn  cp 


(2.4) 


where    m    refers    tc   the   mixed    layer,    s    to    the    surface    value, 


is    the    heig! 


-F     •*-  \ 


he    mixed    layer,    Z    is    the    height,     &<j>    is    the 


jump  in  <j>  across  the  inversion  zone,  r  is  the  gradient 
above  the  inversion,  and  r.  is  the  gradient  in  the  irixsd 
layer    (assumed   tc    be   zero) . 

2-      Dynamics   within   the    mixed   lay er 

The  model  predicts  the  time  evolution  of  the  mixed 
layer  values  of  g  ,  9  ,  and  h  (see  Egn.  2.3).  At  initiali- 
zation, the  lapse  rates  -for  these  quantities  are  held 
constant  while  larger,  synoptic  scale  observations  provide 
the  input  values  of  the  boundary  layer  wind  for  the  entire 
forecast   period. 

Surface  fluxes  of  sensible  heat  and  water  vapor  are 
deteriined  through  an  application  of  the  bulk  aerodynamic 
formula  (Davidson  et  al.  1978).  These  fluxes  are  driven  by 
the    bulk   air/sea    differences   and    are    expressed   as 


{wtQ'Jn  =  <=2    le.  -  e   )U 


e     0 


8  0 


m 


(2.5) 


("X'o   "   ce    (<>o   "  Vu 


where    CQ   is    the    ten    ireter,    stability    dependent,    drag    coeffi- 
cient,     u*    is   the   friction    vel 
respective    sea-surface    values. 


cient,      u*    is   the    friction    velocity,         and    9      and    q      are    the 
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rar ameterizaticn  for  the  entrainment  rate  is  that 
given  by  Stage  and  Eusinger  (1981)  »  This  parameterization 
allows  closure  of  the  system  of  equations  and  determines  the 
time  evolution  of  the  inversion  height.  It  also  brings  the 
turbulent  kinetic  energy  (TKE)  budget  into  the  formulation 
since  entrainment  of  heat  flux  across  the  stable  inversion 
requires  TKE,  that  is,  work  is  done  by  this  process.  Here 
the  assumption  used  fcr  closure  states  that  the  dissipation 
rate  cf  TKE  is  a  fixed  fraction  (1-A)  cf  the  production 
rate,  where  A  is  the  entrainment  coefficient,  set  equal  to 
0.2.  Within  the  TKE  equation  are  terms  which  provide  for 
production  and  dissipation,  therefore  it  is  assumed  that 
both  occur  within  the  mixed  layer,  rather  than  only  where 
there      is   negative      buoyancy      flux.  These    are      associated 

through    the      use   of    the      turbulent    scaling    velocity      for   the 


mixed    layer. 


■*  * 


and   the      mixed   layer    depth    (h)         for    the 


scaling    length. 

3  •      Bad  dative   flux   cons  i derations 

fchen  considering  radiation  flux  and  its  effects  on 
the  mixed  layer  there  are  several  points  which  must  be 
understood.  Distirction    must      be      made    between      long-wave 

fluxes    for    clear   or    cloudy    cases,    and    short-wave    fluxes. 

a.      Long-wave  radiative   flux 

Onder  cloudless  sky  conditions,  the  net  long- 
wave radiative  flux  is  determined  using  the  water  vapor  and 
temperature  profiles.  The  flux  is  calculated  by  inte- 
grating, from  the  surface  (Z=0)  to  the  top  of  the  mixed 
layer  (Z  =  h)  ,  the  flux  emissivity  profile  where  the  emis- 
sivity  is  determined  from  the  Stef an-Boltzmann  law  (Fleagle 
and  Eusinger  1980) .  The  aforementioned  modification  has 
also    beer,   applied. 
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The  presence  of  clcuds  necessitates  the  calcula- 
tion of  the  net  radiative  flux  at  both  the  cloud  tcp  and  the 
cloud    base.  These    calculations      are    also      made    using      the 

Stef an-Ecltzmann   law.  When   the      cloud   top      flux   is      being 

determined,  the  input  temperatures  are  the  cloud  tcp  temper- 
ature and  the  effective  radiative  sky  temperature.  Fcr  the 
cloud  base  the  input  temperatures  are  the  temperatures  at 
the  cloud  base  and  at  the  sea  surface.  It  is  also  important 
to  note  that,  in  this  case,  the  flux  divergence  between  the 
sea    surface   and    the    cloud   base    is   neglected. 

b.      Short-wave   radiative    flux 

The  «ABL  model  calculations  of  short-wave  radia- 
tion take  into  consideration  both  direct  and  diffuse  compo- 
nents. Eoth  scattering  and  absorption  contribute  to 
short-wave  extinction.  The  scattering  ccntributicn  tc  the 
short-wave  radiative  flux  is  determined  through  use  of  the 
delta-Eddington  methcd,  where  the  mixed  layer  scattering  is 
due  tc  atmospheric  particles,  cloud  droplets  in  the  cloudy 
case  and  sea-salt  aerosols  in  the  clear  sky  case. 
absorption  is  primarily  the  result  of  water  vapor  and  is 
taken  into  account  when  the  incident  flux  at  the  top  of  the 
mixed   layer   is    determined. 

C.       HCEEI    AFPLICATIOl 

All    model  runs      for   this   thesis   were   done      such    that   the 
situation   approximated     to      a   real    forecast   run.  That   is, 

the  model  was  initially  set  tc  run  with  a  single  sounding 
and  an  estimate  of  -C.0005  m/s  for  the  subsidence.  Further 
runs  were  made  as  each  new  sounding  could  be  added.  This 
presents  a  realistic  representation  of  what  one  could  expect 
from    the    mcdel    for    the    Vandenberg    AFB   area. 
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Preparation    for    running    the      model   requires   that    sever?.! 
points   or   i^ems    be    considered.      These   consist   of    identifying 

the  synoptic  pattern  and  initialization  parameters;  will  be 
discussed   below. 

1  .      Identify   S vnc^tic   Pattern 

Eefcre  checking  on  input  data  for  the  model  the 
overall  syncptic  pattern  must  be  checked-  To  do  this  one 
utilizes  all  available  meteorological  data  (e.g,  surface  and 
upper  air  synoptic  analyses,  as  well  as  satellite,  rawin- 
sonde  and  individual  surface  observations).  The  desired 
synoptic  situation  fcr  optimum  model  predictability  is  to 
have  the  area  cf  interest  dominated  by  a  high  pressure 
syste  ir . 

2 •      Identify    Initializa tion   Parameters 

The  next  step  in  preparing  to  make  the  model  runs  is 
to   determine      the   input    parameters.  Potential   temperature 

and  moisture  values  fcr  the  mixed  layer,  height  and  strength 
cf  the  jump,  surface  wind,  sea-surface  temperature,  vertical 
velocity  and  the  lapse  rate  above  the  jump  are  the  necessary 
parameters  and,  of  these,  a  good  determination  of  the  jump 
value    across  the   inversion    is   critical. 

Except  for  vertical  velocities,  parameters  are 
determined  from  the  vertical  soundings,  and  can  be  done 
manually  cr  through  the  use  of  a  digitizing  program.  k 
sample  cf  a  digitized  rawinscnde  appears  in  Fig.  2.2. 
Extreme  care  must  te  used  sc  that  the  values  determined 
represent   the   true    scunding    or   serious   error    can    result. 

The  vertical  velocity  or,  in  our  case,  the  subsi- 
dence rate  for  the  forecast  start  time  must  be  determined. 
There  are  several  means  for  calculating  vertical  velocity 
but  ncne,  to  date,  provide  the  accuracy  desired.  with  this 
in   mind,    the  method    which   has    proved   to    be    most    accurate    was 
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used.         This   method      is   the    Q-irethod    which      is  an    expression 
for    vertical  velccitj   derived      from   the    vertical    integration 


cf   the    cccisture    budget   equation    and    is   given    by: 


A  (q  +  1)   (dh/dt)     -     h(d(q  +  l)/dt)      +     (w'q1) 


w    = 


0 


A(q+1)     +     3(h/2) 


(2.6) 


where  g  and  1  are  vapor  and  liquid  water  contents,  ana  8, 
their  vertical  gradient  (Gleason,  1982).  The  computer 
program  used  to  calculate  the  vertical  velocity  uses  current 
data  as  sell  as  data  within  the  previous  24  hours.  The  use 
cf  previcus  rawinsonde  data  is  allowed  since  -che  Q-method 
equation  includes  total  change  with  raspect  to  time  (d  /dt) 
for      th =   inversion      height      and      the   total      moisture.  The 

vertical  velocities      cutput    frcm      the    program     are    for      each 
data    period   input.  These    values   are  then      averaged   to   get 
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Figure  2.2    Sanple  of  digitized  rawinsonde 
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the  rate  fcr  the  ertire  pericd  with  the  final  sutsicer.ee 
rate  calculated  teinc  the  one  used  in  the  model*  The  only 
difficulty  encountered  with  this  method  was  when  the  subsi- 
dence inversion  height  was  increasing  at  a  rata  sufficient 
to  yield  a  positive  average.  When  this  occurred  the  mcdel 
input  rate  was  set  tc  zero.  Table  I  gives  the  subsidence 
rates    (i/s)    used   for   the   model   runs. 


TABLE    I 
Subsidence  Hates   for   Model   Runs 


DATE 

GMT 

SDESIDNECE 

DATE 

July 

25 

00 
12 

-0.C05 

-0.003948 

August 

26 

00 
12 

-0.003359 
-0.0008265 

27 

00 
12 

-0.01489 
-0.CO1553 

28 

00 

-0.002856 

September 

December 

13 

00 
12 

-0.C05 
-0.005813 

14 

00 
12 

-0.C02856 
0.0 

April 

13 

00 
12 

-0.C05 

-O.C02019 

14 

00 
12 

O.C 

-0.000365 1 

15 

00 

12 

O.C 
O.C 

16 

00 
12 

-0.C02208 
-0.002988 

GMT  SUBSIDENCE 

29  00  -0.005 

12  -0.0C08393 

30  00         0.00 

12  -0.0002389 

31  00  -0.0004126 
12         0.0 

1  00         0.0 

2  00  -0.005934 
12  -0.001676 

3  00  -0.000082 
12         0.0 

4  00  -0.0006649 
12  -0.0002713 

5  00  -0.0005856 
12  -0.001524 

6  00         0.0 
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III.     SYNOPTIC    PATTEBN    DISCISSION 

Fcur  data  sets  vera  used  for  evaluation  of  model 
validity.  These  sets  include  a  summer  period  (25  through  27 
July  1977),  a  winter  period  (13  and  14  December  1977) r  a 
sprirg  period  (13  through  17  April  1981),  and  a  fall  period 
(29  Aug  through  5  September  1981).  Large  scale  analyses 
were  for  as  close  tc  1200  Greenwich  Mean  Tims  (GMT)  as 
possible.  In  most  cases  the  satellite  pic-ures  will  te  at 
17U5  GM1  since  these  pictures  best  display  the  area  of 
interest . 

This  discussion  will  be  of  the  overall  weather  picture 
for  the  eastern  North  Pacific  Ocean  and  the  western  United 
States.  Included  will  be  conditions  over  the  area  as  a 
whole  and  unique  or  differing  conditions  at  Vandenberg  AFB 
itself.  Satellite  pictures,  surface  weather  analyses,  and 
upper  air  analyses  (5C0  millibars)  will  be  used  tc  determine 
the    synoptic  situation  of   each    period. 

1.       SUMMER    CASE,    PESICD    I    (25-27      JULY.    1977) 

The  surface  and  500  millibar  (mb)  synoptic  patterns  for 
this  period  appear  in  Figs.  3.2  through  3.4.  During  this 
pericd  the  coastal  region  was  dominated  by  a  high  pressure 
systeir  ever  the  North  Pacific  Ocean  with  the  thermal 
troughing  extending  from  Mexico  northward  into  the  high 
desert  regicn  of  Arizona  and  southeastern  California.  Such 
a  synoptic  pattern  usually  produces  weak  onshore  flew  along 
the      southern   California      coastal   regions.  The    upper      air 

analyses       (500    mb)  showed      no      strong    westerly       flew      over 

California    during   the   entire   period. 


2a 


During  this  period  the  Pacific  high  pressure  system 
moved      northeastward      and      on      27  July,         weakened      as      a 

stationary  front  developed  off  the  northwest  coast.  Also  of 
importance  en  27  July  was  the  development  of  a  low  pressure 
center  west  of  Los  Angeles,  California  Even  with  this 
development  the  predominate  flow  remained  onshore  in  the 
?anderbeig    AFB    region. 

The  rawinsonde  data  revealed  a  subsidence  type  inversion 
throughout  the  period,  see  Fig.  3.1.  There  was  also  a 
diurnal  variation  in  the  inversion  base.  Each  day  the  base 
ranged  from  the  surface  at  1200  GMT  to  a  height  above  650 
feet  at  COOC  GMT.  An  example  of  this  diurnal  change  also 
appears    in    Fig.     3. 1 . 

The  satellite  pictures  for  this  period  appear  in  Figs. 
3.5  through  3.7.  Orly  patchy  areas  of  lew  clouds  along  the 
Califcrria  coast  appear  on  the  satellite  pictures.  These 
areas  were  in  the  San  Francisco  bay  area  and  the  coastline 
south  of  Lcs  Angeles.  The  only  ether  area  of  low  clouds  i»as 
two  tc  three  degrees  longitude  west  of  Vandenberg  AFB. 
There  was  no  significant  changes  to  the  coastal  locations  of 
the  lew  clcuds  while  the  area  over  the  water  had  some  dissi- 
pation and  was  even  further  eff shore.  The  Vandenberg  AFB 
vicinity  remained  under  clear  skies  during  the  entire 
period. 
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Figure  3-1   _ Example  of  Summer  Case  rawinsonde  and 
diurnal  cnange  in  inversion. 
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Figure  3. 2 


Surface    (a)    and   500    mb    (b)    analyses 
1200   G3T    25      July    1977. 
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Figure  3.3    Surface  (a)  and  500  mb  (b)  analyses 
1200  GMT  26   July  1977. 
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Surface  (a)  and  500  ab  (b)  analyses 
12C0  GHT  27   July  1977. 
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Figure  3.5    Satellite  picture  for  25   July  1977. 
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figure  3.6   Satellite  picture  for  26  July  1977. 
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figure  3-7    Satellite  picture  for  27   July  1977. 
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E.       ilN'IEB    CASE,    PEB2CD    II     (13-14    DECEMBER    1977) 

The  general  surface  and  500  mb  analyses  appear  in 
figures  3.S  and  3.10  A  difference  between  this  syncptic 
patterr  and  that  found  in  the  other  periods  is  in  this 
pattern  th  =  desired  onshore  flow  exists  but  it  does  not 
arise  frcm  the  strong  Pacific  high  and  thermal  troughing  tut 
from  a  high  pressure  cell  that  had  separated  frcm  the  main 
Facific  cell  and  moved  in  over  central  California  The 
strongest  westerly  flow  aloft  was  found  to  the  north  of 
Califcrria  with  the  Polar  jet  ircving  across  the  Washington  - 
Cregcc      region.  Westerly      flow   over      Vandenberg      AFB      was 

stronger   than  in   other   periods. 

The  primary  aspect  is  that  on  13  December  a  high  pres- 
sure center  moved  over  central  California  with  the  large 
high  pressure  system  still  maintained  over  the  North  Pacific 
Ocean.  Weak  thermal  troughing  was  also  extending  northward 
into  California  The  overall  result  of  these  surface 
patterns  *as  weak  onshcre  flew  in  the  Vandenberg  AFB  and  San 
Nicolas  Island  regions.  The  enly  other  significant  feature 
in  the  13  December  analysis  was  the  frontal  system  mcving 
into    the    Pacific    Northwest. 

Cn  14  December  the  high  pressure  center  dissipated 
leaving  the  Pacific  high  pressure  system  and  ridging  acrcss 
central  California.  Weak  thermal  troughing  continued  in  the 
southern  California  area.  Offshore  flow  also  continued  to 
dominate  the  southern  Califcrria  coastline.  The  frenfal 
system  frcm  the  13th  moved  eastward  and  a  trailing  warm 
front    moved  to    northern   Califcrria. 

Fig.  3.8  is  an  example  of  the  vertical  soundings  for 
this    period.  The    rawinsonde      data    showed      subsidence   type 

inversicrs  on  both  days.  Most  important  is  that  on  13 
December  at  0000  GMT  the  inversion  was  very  weak  and  not 
well    structured.  The    inversion    strengthened      and    remained 
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aloft    until    1200   GMT   en    14      December    when   it   was    observed   as 

a      surface    tased      inversion.         At      this      rime    the      inversion 
reached    the    maximum    strength    for    the   period. 

The  satellite  pictures  for  this  period  appear  in  rigs. 
3.11  and  3.12.  Clear  skies  in  the  Vandenberg  AFB  area  with 
a  frontal  system  movirg  intc  the  northern  California  regions 
appear  in  the  satellite  picture  for  13  December.  There  is 
very  little  indication  of  low  clouds  along  the  southern 
California  coastline.  By  14  December  the  front  is  well  into 
northern  California  with  increased  cloudiness  across  central 
California .  There  has  also  been  an  increase  in  the  amount 
of  lew  clouds  along  the  southern  California  coast.  Patchy 
stratus  fellows  the  coast  south  of  Los  Angeles  and  a  large 
area    now   sits   off  the   Vandenterg   AFB    coast. 
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Figure  3.11    Satellite  picture  for  13  December  1977. 
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Figure  3,12    Satellite  picture  for  14  December  1977. 
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C.       SPRING    CASE,    PERICD    III     (13-16    APRIL    1981) 

The  syncptic  patterns  for  this  period  appear  in  Figs. 
3.14  through  3.17.  The  dominating  feature  was  a  high  era- 
sure system  extending  northeastward  from  the  North  Pacific 
Ccean      regions      across      Washington      into      Canada.  Therial 

trougfcinc  dominated  the  high  desert  regions  of  California 
and  extended  into  northern  California.  The  predominate 
southern     California      coastal      flew,  generated      by      these 

patterns,  was  onshore.  500  mb  analyses  shew  light  and  vari- 
able winds  with  tendencies  toward  formation  of  a  lew  pres- 
sur a  center  through  the  period.  The  strongest  westerly  flow 
remained    well   north    of   the    Vandenberg    AFB    area. 

Cn  13  April  the  Pacific  high  was  strong  but  showing  a 
tendency  to  weaken  as  the  longwave  pattern  extended  into 
Canada.  As  the  longwave  high  pattern  shifted  to  the  central 
Dnited  States  the  Pacific  high  weakened.  The  thermal  trough 
in  southern  California  continued  to  dominate  but  alsc  had 
some    weakening. 

Ey  15  April  the  Pacific  high  system  had  begun  to 
strengthen  with  some  ridging  across  central  California.  The 
high  continued  intensification  through  the  remainder  of  the 
period.  Cn      16      April     the        ridging      extended      into      the 

Washington  area  with  the  Pacific  high  pressure  center  begin- 
ning to  ireve  northeastward.  Throughout  this  portion  of  the 
period  the  thermal  trough  showed  very  little  variation.  The 
only  ether  key  occurrence  from  the  15th  to  the  16th  was  the 
movement  of  a  cold  front  across  the  Pacific  Northwest  into 
Nevada. 

Fig.  3.13  is  an  example  of  the  rawinsonde  soundings  for 
this  period.  Rawinsonde  retorts  during  this  period  had 
subsidence    inversions      with    much      daily    variation.  At   the 

beginning  of  the  period  the  inversion  was  weak  and  aloft. 
Diurnal      variations      intensified      the         inversion      but      also 
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brought  it  to  the  surface.  This  pattern  continued  through 
the  14th  but  on  the  15th  the  inversion  intensified  at  3G00 
GMT  and  remained  alcft  until  1200  GMT  on  the  16th.  During 
the  latter  half  cf  the  period  there  was  some  fluctuation  in 
the  intensity  but  over  all  the  inversion  remained  strong. 

These  satellite  pictures  appear  in  Figs.  3.18  through 
3.21.  On  13  April  the  satellite  showed  Vandenberg  AFE  under 
clear  skies  while  the  coastal  region  south  of  there  was 
dominated  by  low  clouds.  Ey  14  April  the  area  cf  coverage 
ty  low  clouds  increased  and  Vandenberg  AF3  was  within  that 
area.  low  clouds  continued  to  dominate  the  Vandenberg  AFB 
area  with  seme  dissipation  beginning  on  15  April.  16  Apr 
showed  criy  thin  coverage  of  low  clouds  over  the  Vandenberg 
AFB  area  with  the  denser  coverage  south  of  the  area. 


3800r 


2500' 

2000- 

1500 

1000- 

500' 
100- 


THETTfl 


850- 


■1000- 


SPGC.    HUM. 

(g-'Wg) 

5 

10 

15 

20 

25 

POT.     TETMP. 

to 

0 

10 

20 

30 

40 

Figure  3.13    Example  of  rawinsonde  for  Spring  Case 
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Fiqure  3.14    Surface  (a)  and  500  mb  (b)  analyses 
1200  GST  13  April  1977. 
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(a) 
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figure  3,15    Surface  (a)  and  500  mb    (b)  analyses 
12C0  GMT  14  April  1977. 
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Jigure  3.16 


Surface  (a)  and  500  mb  (b)  analyses 
1200  GHT  15  April  1977. 


U3 


(a) 


r        \      ^Vi 


13-Sr5    .' 


(b) 


Figure   3.17 


Surface    (a)    and   500    mb    (b)    analyses 
1200   GHT    16    April    1977. 
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Figure  3.18    Satellite  picture  for  13  April  1981. 
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1745  14APS1  36A-2   00243  16021  US2 


figure  3.19    Satellite  picture  for  14  April  1981. 
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1745  1 5 APS 1  36A-2   00243  16002  WB2 


Iigur€  3.20    Satellite  picture  for  15  April  1981. 
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1745  16APS1  36A-2   00251  15992  U82 


Figure  3.21    Satellite  picture  for  16  April  1981. 
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B.       JILL    CASE,     PERIOD   IV    (29    AOGtJST    -    5    SEPTEMBER    1981) 

Ihe  surface  {a)  and  50  0  m fc  (b)  analyses  for  Period  17 
appear  in  Figs.  3.2  3  through  3.30.  The  predominate  surface 
feature  during  this  extire  period  was  the  presence  of  a  high 
pressure  system  in  the  Eastern  portions  of  the  North  Pacific 
Cceari.  A   thermal      low     and      trough   system      dominated      the 

southern  California  region.  The  effect  of  these  two  surface 
features  was  to  produce  weak  onshore  flow  in  the  area  of 
concern.  Conclusions  drawn  from  the  500  mb  analyses  show 
that,  except  for  the  last  day  of  the  period,  the  maximum 
westerly      flow    remained     north   of      California.  Only    en      5 

September  did  the  Pclar  jet  extend  toward  the  regions  of 
interest. 

The  surface  analysis  for  29  August  (Fig  3.23)  shows  the 
western  Pacific  Ocean  dominated  by  a  high  pressure  system 
with  a  lew  pressure  center  located  off  the  Washington  coast. 
A  thermal  low  extended  northwestward  from  Mexico  into 
Calif crria        with        trough! ng  extending        into        northern 

California.  The  southern  end  of  a  cold  front  had  also  moved 
into    the   northern  California   area. 

Ey  30  August  the  cold  front  had  moved  into  southern 
California  and  the  thermal  low  had  changed  its  orientation. 
The  Pacific  high  pressure  system  showed  weak  intensification 
with  irevement  toward  the  northeast  and  had  pushed  inland  by 
the    3  1st.  The   31st      also    shewed   that      the    cold      front    was 

beginning  to  leave  California  but  with  another  frontal 
system  iroving  out  cf  the  Gulf  of  Alaska.  Figures  3.24 
through    3.25  depict    these   occurrences. 

The  1  September  analysis  revealed  a  weakening  of  the 
high  pressure  and  a  cold  front  approaching  California.  A 
thermal  trough  was  well  established  over  California  tut  very 
diffuse    en    2  September    as   the    cold    front    moved   through 
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southern  Calif orria.  The  high  pressure  system  agair.  inten- 
sified and  showed  no  marked  weakening  throughout  the 
remainder  cf  the  period.  The  only  significant  event  after  2 
September  was  the  entrance  of  a  rapidly  moving  ccld  front 
into  California.  This  scenario  is  shown  in  Figs  3.26 
through    3.30. 

A  sairple  of  the  rawinscndes  for  this  period  appears  in 
Fig-  3.22.  The  rawinscndes  for  this  period  have  a  very  good 
subsidence    type      inversion    throughout    the    period.  In   each 

rawinscr.de  release  the  inversion  was  found  above  1000  feex 
with  decreasing  relative  humidity  through  the  inversion 
zone. 

The  extent  and  development  of  the  low  clouds  appear  in 
the  pictures  in  Figs.  3.31  and  3.38.  The  satellite  pictures 
revealed  a  buildup  cf  coastal  low  clouds  (stratus  and/or 
fog)  during  the  period.  Beginning  on  29  August  the  clouds 
tended  to  te  thin  with  seme  dissipation  during  the  day. 
Further  intc  the  period  the  stratus  thickens  and  by  the  end 
of  the  period  there  is  a  dense  area  of  low  clouds  off  the 
central   and    southern   California   coastal    regions. 
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Figure  3.22        Bawinsonde   example   for   Fall   Case. 
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Jigure   3.23 
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Scrface  (a)  and  500  ib  (b)  analyses 
12C0  GHT  29  august  1977. 
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(a) 
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Eigure  3.24 


Surface  (a)  and  500  mb  (b)  analyses 
1200  GHT  30  august  1977. 
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Figure  3,25 


Surface  (a)  and  500  mb  (b)  analyses 
1200  GHT  31  August  1977. 
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liqure  3-26    Surface  (a)  and  500  mb  (b)  analyses 
1200  GHT  1  September  1977. 
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Piqure   3.27        Surface     (a)    and   500    mb    (b)    analyses 
1200   GHT   2    September    1977. 
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Figure   3.28 


Surface     (a)    and   500   mb    (b)    analyses 
1200   GHT   3    September    1977. 
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Fiqure  3.29    Surface  (a)  and  500  ab  (b)  analyses 
1200  GHT  4  September  1977. 
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Figure  3.30    Surface  (a)  and  500  mb  (b)  analyses 
1200  GHT  5  September  1977. 
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500  29AU81  17A-2 


figure  3.31    Satellite  picture  for  29  August  1981. 
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Figure  3.32   Satellite  picture  for  30  August  1981. 
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Figure  3.33    Satellite  picture  for  31  August 
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1745  0ISE81  36A-4       00532  19161"  UC2 


Figure   3.34        Satellite    picture   for    1    September    1981. 
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17-ro  02SE81   3bA-2       002.4  16002  MB2 


Figure  3.35   Satellite  picture  for  2  September  1981. 
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Figure  3.36    Satellite  picture  for  3  September  1981. 
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Figure  3.37    Satellite  picture  for  4  September  1981. 
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Figure  3.38    Satellite  picture  for  5  September  1981 
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IV.  RESULTS 

Model  forecasts  were  made  for  every  30  minutes  for  24 
hours  following  the  initialization  time.  The  forecasts 
presented  in  the  following  tables  show  only  the  12  and  24 
hour  values  since  they  correspond  to  the  verification  times. 
The  structure  of  the  tables  are  such  that  the  first  fore- 
cast, indicated  by  column  heading  FCST1,  was  based  on  the 
first  observation  for  initialization  and  within  that  cclumr 
are  the  forecasts  generated.  For  example,  if  the  first 
observation  was  for  1/0000  GMT  then  the  forecasts  from  that 
observation  would  be  for  1/1200  GMT  and  2/0000  GMT.  These 
forecasts  would  be  iv.  column  FCST1  .  The  second  forecast, 
column  heading  FCST2,  was  initialized  with  the  second  obser- 
vation (e.g.,  1/1200  GMT)  and  within  that  column  ar?  the 
forecasts  generated  (2/0000  GMT  and  2/1200  GMT).  Each  table 
continues  ir.  this  manner  until  all  forecasts  for  a  period 
have  beer  given. 

Boot  mean  sguare  (RMS)  scores  are  given  following  each- 
forecast  except  the  last.  No  EMS  for  the  last  forecast  were 
given  as  there  were  ec  observations  for  comparison.  Also 
note  that  the  12  hour  RMS  value  simply  represents  the 
difference  between  observed  and  forecasted  values. 

The  values  given  in  the  tables  ar  =  for  the  parameters: 
inversion  height,  lifting  condensation  level,  mixed  layer 
potential  temperature  and  mixed  layer  specific  humidity. 
Units  for  these  are  ireters  for  the  inversion  height  and  the 
lifting  condensation  level,  degrees  Centigrade  for  potential 
temperature,  and  g/kc  for  specific  humidity.  Also  note  that. 
subsidences,  used  in  discussing  these  results,  are  given  in 
Table  I. 
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3.       2ISCCS5T0N    OF    EOEEL    RESULTS    FOR    PERIOD    I     (SUMMER) 


This    period    was    dominated 


clear    sfcxes. 


rpV   - 


lifting   condensation    level    values    (see   Table    III)     only 


fore- 


casted     clcuds       during      FCST2.  The      forecasted      inversion 

heights,  Table  II,  reveal  that  the  model  was  either  very 
close  or  quite  in  error  relative  to  the  observed  height. 
Tables  IV  and  V  reveal  that,  throughout  the  period,  there 
were    alsc    problems    associated    with   potential  temperature    and 


TABLE    II 
July   Inversion   Results 


DATE    GET    CBS 
25  00    356 


FCST1       RMS       FCST2       RMS       FCST3       RMS 


12 

157 

195 

38.0 

- 

- 

- 

- 

26 

00 

277 

108 

173.2 

92 

185.0 

- 

- 

12 

199 

- 

- 

112 

204.4 

165 

3  4.0 

27 

00 

206 

— 

— 

— 

— 

100 

111.3 

DATE 

G^T 

C3S- 

FCST4 

RMS 

FCST5 

RMS 

FCST6 

RMS 

FCSI7 

2  7 

CO 

206 

172 

34.0 

- 

- 

- 

- 

12 

120 

146 

42.8 

154 

34.0 

- 

- 

- 

28 

CC 

199 

- 

- 

116 

89.7 

71 

128.0 

- 

12 

- 

- 

- 

- 

- 

42 

- 

1  10 

2S 

CO 

- 

- 

- 

- 

- 

- 

62 

TAELE    III 
July    Lifting    Condensation    Levels 


ATE 

GMT 

OBS 

FCST1 

ECST2 

FCST3 

FCST 

25 

00 

633 

- 

- 

- 

— 

12 

146 

2  99 

- 

- 

- 

26 

00 

1119 

199 

42 

- 

- 

12 

217 

- 

0 

519 

- 

27 

00 

206 

- 

- 

309 

185 

12 

452 

- 

— 

— 

168 

ATE 

GMT 

OBS 

FCST5 

FCST6 

FC3T7 

27 

12 

452 

605 

- 

- 

28 

00 

1261 

410 

184 

- 

12 

- 

- 

176 

586 

29 

00 

- 

- 

- 

366 
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TABLE    I? 
July   Potential   Temperatures 

DATE    GMT       OBS      FCST1  FCST2  FCST3         FCST4 


25 

0  0 

16.2 

- 

- 

- 

- 

12 

11.  6 

17.4 

- 

- 

- 

26 

00 

17.  4 

13.2 

11.6 

- 

- 

12 

12.6 

- 

8.5 

17.8 

- 

27 

00 

20.9 

- 

- 

18.4 

15.  2 

12 

14.  1 

— 

— 

— 

16.7 

DATE 

Gai 

GB5 

FCST5 

F  C  5  T  6 

FCSI7 

27 

12 

14.  1 

19.3 

- 

- 

28 

00 

19.5 

19.6 

17.  1 

- 

12 

- 

- 

18.3 

18.9 

29 

00 

- 

- 

- 

19.3 

TABLE    V 
July  Hoisture   Values 

DATE    GMT  OBS  FCST1  ECST2          FCST3         FCST4 

25  00  8.2  - 

12  7.9  6.3  - 

26  00  6.8  5.4  3.3 

12  8.1  -  7.2               9.7 

27  00  9.3  -  -               11.3               9.8 
12  7.9  -  10.8 

DATE    GMT  OBS  FCST5  ECST6          FC5T7 

27  12  7.9  10.6 

28  00  7.3  11.6  11.1 

12  -  -  12.0            10.1 

29  00  -  -  11.6 


specific   humidity    forecasts.      These    problems    will    be   covered 
during    the    following    discussions. 

Examination  of  the  inversicn  height  results  reveals  seme 
good  and  seme  poor  height  fcrecasts  at  the  12  hour  pcint  bur 
mostly  bad  forecasts  at  24  hours.  Those  forecasts  made  with 
average  sutsidence  rates  (-0.003  to  -0.005  m/s)  had  fair  to 
good  results  at  12  hcurs  but  had  poor  results  at  24  hours. 
The  bes-  height  forecasts  (Table  II  column  FCSTU)  were  with 
the  weakest  subsidence  rate  while  the  second  best  (The  fifth 
forecast)  vas  with  the  strongest  subsidence  rate.  Also  of 
interest      is  that      fcrecasts    with      the   strongest      subsidence 
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also  had  the  lifting  condensation  level  well  abcve  th;- 
inversicr.  Another  important  feature  in  Table  II  is  that, 
except  for  the  forecast  with  the  weakest  subsidence,  the 
best  forecasts  were  for  the  12  hour  periods  initialized  at 
1200    6H1    (a   morning    period    for   the   California   region). 

Potential  temperature,  specific  humidity,  radiation 
(with  and  without  clouds)  and  entrainment  will  be  examined 
in  explaining  what  was  occurring  with  the  predicted  ir.ver- 
sion      heights    for      this      period.  Each    parameter      will      be 

described  separately  as  to  its  possioie  impact  on  the  inver- 
sion   forecast. 

1)  potential  temperature  -  Under  clear  skies  one 
would  normally  expect  nightly  cooling  ever  land. 
The  observed  values  shown  in  Table  IV  reveal  that 
this  occurred  yet  the  forecasts  do  net  reflect  this 
fact.  Nocturnal  cooling  strengthens  the  jump  and 
when  coupled  with  subsidence  lowers  the  inversion  as 
seen  in  Fig.  4.1.  All  forecasts  were  lowering  the 
inversion  due  to  good  subsidence  rates  for  initiali- 
zation. Daytime  heating  acts  against  the  subsi- 
dence, weakens  the  jump,  and  results  in  lifting  of 
the  inversion.  This  diurnal  event  was  not  reflected 
in    the    forecasts. 

2)  specific  humidity  -  After  the  first  two  fore- 
casts, specific  humidity  forecast  values  were  only 
fair.  There  were  nc  forecast  values  which  were  good, 
since  the  errors  ranged  from  just  under  1  g/kg  up  to 
3.3  g/kg.  The  primary  trend  in  the  forecasts  was  to 
increase  specific  humidity.  This  prediction  of  too 
large  a  specific  humidity  caused  the  model  to 
predict  even  lower  lifting  condensation  levels 
which,  in  turn,  enhances  the  poor  inversion 
predictions. 


71 


3)  radiation  -  This  factor  must,  be  considers.?  under 
conditions   of  clear    as    well    as    cloudy   skies, 

a)  cloudy  -  In  the  second  forecast  the  model 
predicted  clouds  to  for  a  and  extend  to  the 
surface.  The  clouds  radiate  long  wave  energy  to 
space  so  the  potential  temperatures  were  forecast 
too  low.  Inversion  predictions,  while  decreasing 
then  increasing  during  this  forecast,  did  not 
reflect  the  diurnal  change  predominate  with  clear 
skies . 

h)  clear  -  Radiation  computations  here  could 
have  been  incorrectly  influenced  by  the  large 
specific  humidity  forecasts.  This  accompanied 
with  good  subsidence  rates  caused  all  fcrecasts 
to   lower    the    inversion. 

4)  entrainment  rates  -  these  influence  the  time 
rate  of  chance  of  specific  humidity  and  potential 
temperature  as  well  as  inversion  height.  In  several 
of  the  forecasts  the  entrainment  rate  was  zero  so  it 
was  internally  set  tc  a  positive,  yet  minimal, 
value.  This  could  have  affected  the  ether  fore- 
casted values  but  it  is  believed  the  effect  was  less 
than  others. 

A   definite    diurnal   pattern   can   be    seen    (see    Fig.    4.1)    in 
the      cfcserved      inversion      heights.  Generally      there      were 

increased  heights  in  the  transition  from  morning  to  evening 
and    decreases      fcr   the      opposite   transition-  This   pattern 

coincides  with  the  diurnal  pattern  of  the  potential  tempera- 
tures which  reflect  radiative  daytime  heating  and  nightly 
cooling.  Surface    winds      also   displayed      the  influence     of 

daily  heating  with  evening  reports  having  the  strongest  wind 
speeds  and  with  steady  northwesterly  directions  (an  onshore 
wind).  Figure  4.1  also  depicts  the  12  and  24  hour  fcrecasts 
and    it    is   easily   seen    that      the    diurnal    pattern    shown    by  the 
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observed  data  is  net  shewn  ry  the  forecasts.  Only  the 
morning  tc  evening  12  hour  forecast  indicate  a  correct 
trend.  However  close  inspection  of  the  forecasts  indicate 
that  at  these  times  the  model  inputs  for  subsidence  were 
strong  thereby  influencing  the  forecast. 

Thsse  effects  all  suggest  problems  with  the  radiation 
calculations  throughout  this  summer  period.  A  reason  for 
this  error  is  the  fact  that  Vandenberg  AFB  is  a  shoreline 
station  and  the  model  is  primarily  designed  for  over  xafer 
use.  Therefore,  land  heating  and  cooling  by  solar  radia- 
tion, even  through  a  stratus  sky,  reguire  adjustments  to 
apply  the  model  at  a  shoreline  station.  Vertical  velocity 
and  entrainment  rates  also  affected  forecasted  MABL  proper- 
ties through  influences  on  the  rate  of  change  of  key  parame- 
ters as  well  as  the  inversion  heights  (see  Eqn.  2.3) . 

B.   DISCISSION  OF  MODEL  RESDLTS  FOR  PERIOD  II  (WINTER) 

Besults  in  Table  VI  generally  suggest  that  during  -his 
season  the  model  fails  to  forecast  variations  in  the  height 
cf  the  inversion.  Tables  VII  through  IX  compare  forecast 
and  observed  lifting  condensation  levels,  potential  tempera- 
tures and  specific  humidity  values. 


TABLE  ¥1 
December  Inversion  Results 

DATE  GBT  CBS   FCST  1   RMS   FCST2   RMS   FCST3   RMS   FCST4 
13 

14 
15 


CO 

164 

- 

— 

- 

- 

- 

— 

— 

12 

171 

45 

126.0 

- 

- 

- 

- 

- 

00 

284 

14 

210.7 

40 

244.0 

- 

— 

- 

12 

297 

- 

- 

11 

265.8 

187 

110.0 

- 

00 

- 

— 

- 

- 

- 

124 

- 

492 

12 

- 

— 

- 

- 

- 

- 

- 

623 
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TABLE  VII 
December  lifting  Condensation  Levels 

ME    GMT   OSS   FCST1    FCST2    FCST3    FCST4 
13    00  1083 

14 

15 


12 

20  5  5 

123 

- 

- 

- 

00 

1056 

22 

936 

- 

- 

12 

174 

- 

338 

511 

- 

00 

- 

- 

- 

255 

0 

12 

- 

- 

- 

- 

0 

TABLE  YIII 
December  Potential  Temperatures 

E2TE  GMT   OBS   FCST1    FCST2    FCST3    FCST4 
13 

14 

15 


00 

17.2 

- 

- 

- 

- 

12 

15.8 

14.6 

- 

- 

- 

00 

15.7 

14.  1 

15.4 

- 

- 

12 

13.0 

- 

15.2 

14.9 

- 

00 

- 

- 

— 

14.8 

8.2 

12 

- 

- 

- 

- 

5.4 

TABLE    IX 
December    Hoisture   Values 

ESTE    GMT       OBS      FCST1  FCST2         FCST3         FCST4 

13 

14 

15 


00 

8.9 

- 

- 

- 

- 

12 

3.6 

9.7 

- 

- 

- 

00 

6.3 

9.9 

6.6 

- 

- 

12 

8.  3 

- 

9.0 

8.0 

- 

00 

- 

- 

- 

9.2 

7.3 

12 

- 

- 

- 

- 

7.0 

Vandenberg  AFB  remained  under  clear  skies  in  this  period 
yet  the  forecasts  were  for  rapid  formation  of  clouds. 
Forecast  specific  huiidity  values  improved  in  the  second 
forecast  but  were  still  not  acceptable  for  the  entire 
period.  The  potential  temperature  forecasts  were  best  at 
the    12   hour    point. 
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The  primary  causes  for  the  poor  performance  of  the  model 
in  this  case  appear  tc  re  vertical  velocity  and  entrainsnent. 
The  radiation  and  cloud  calculations  contributed  to  the 
error,  tut  to  a  lesser  extent.  Entrainment  was  a  key  factor 
during  the  first  three  forecasts  because  null  values  were 
obtained  in  computation,  sc  uinimal  positive  values  wsre 
assigrsd.  The  radiation  segment  calculations  are  strongly 
influenced  by  the  presence  cf  clouds  and  on  the  last  day  the 
model  forecast  clouds.  With  zero  vertical  velocity,  the 
model  forecast  rapid  lifting  of  the  inversion.  Kith  no 
observations  it  is  hard  to  evaluate  those  forecasts  tut  the 
trend  was  for  increasing  heights  although  not  as  rapidly  as 
predicted. 

forecast  and  observed  heights  for  this  period  appear  in 
Figure  4.2.  The  diurnal  variations  noted  in  the  summer  case 
are  net  present  in  this  case.  It  is  believed  that  the  large 
errors  in  the  inversion  height  forecasts  are  due  to  two 
factors.  first,  the  vertical  velocities  with  'which  the  runs 
were  initialized  were  large  which  lead  to  forecasts  of 
lowering  heights  in  all  but  the  final  run.  In  reality,  the 
heights  shewed  a  continual  increase.  The  vertical  velocity 
used  in  the  last  run  was  such  that  rapid  lifting  of  the 
inversion  was  forecast.  Seccnd,  the  surface  winds  which 
stayed  strong  even  through  the  morning  hours  allowed  mere 
mixing  in  the  lewer  layers  thereby  increasing  the  depth  of 
the      nixed    layer.  This      is      especially   noticable      in      the 

13/1200    GMT   surface    wind    (360    degrees   at   2.1    m/s)  . 

At  Vandenberg  AFE  a  small  chain  cf  mountains  exists 
due-north  of  the  base.  With  this  northerly  flew,  there 
would  te  an  increase  in  the  local  turbulence  in  the  lower 
levels.  Therefore,  the  poor  height  forecasts  in  this  case 
is  associated  with  the  land  influence  which  is  net  taken 
into    acccunt  by    the    ircdel. 
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C.       DISCUSSION    OP    MODEL    RESULTS    FOR    PERIOD    III    (SPRIHG) 


Tables  X  through  XIII  list  the  forecasted  values  for  the 
spring  cas*;.  This  period  was  cloud  free  the  firs-  day  and 
had  lew  clouds  the  second.  The  forecasts  did  not  include 
cloud  formation  until  the  third  day.  Even  though  the  cicud 
forecasts  were  in  error  it  dees  not  appear  that  the  cloud 
calculations  had  an 3  real  effect  on  the  poor  inversion 
height    predictions. 

The      vertical      velocities      and     entrainment      rates      were 


EAIE    GMT    OBS 
13         00    141 


TABLE    X 
April   Inversion   Results 

FCST1       RMS       FCST2       RMS       ECST3 


RMS 


12 

136 

31 

105.0 

- 

- 

- 

- 

14 

00 

289 

9 

299.0 

75 

214.0 

- 

- 

12 

150 

- 

- 

40 

240.6 

293 

143.0 

15 

00 

469 

— 

- 

— 

- 

298 

222.9 

DATE 

GMT 

OBS 

FCST4 

RMS 

FCST5 

RMS 

FCST6 

RMS 

15 

CO 

469 

140 

329.0 

- 

- 

- 

- 

12 

509 

128 

503.4 

624 

115.0 

- 

- 

16 

00 

421 

- 

- 

835 

429.7 

6  84 

263.0 

, 

12 

291 

- 

- 

- 

- 

867 

663.2 

EATE 

GMT 

OBS 

FCST7 

RMS 

FCST8 

16 

12 

291 

340 

49.0 

- 

17 

00 

- 

440 

- 

323 

12 

- 

- 

- 

335 

different   than      those   for   other    period.  Entrainment    rates 

were  zerc  fcr  the  first  and  third  forecasts,  and  the  first 
three  hours  of  the  fifth  forecast  and  were  assigned  a 
liniiral  positive  value  of  .001.  Hence,  the  inversion  height 
was    subsidence      controlled    for      these    times.  The    moisture 

forecasts  were  good  since  values  agree  to  within  2  g/kg  of 
the  observed  value  and  in  some  cases  even  less  than  1  g/kg. 
Each  of  these  forecasts  were  initialized  with  the  evening 
soundirg  in  which  lifting  condensation  levels  indicated  no 
clouds. 
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TABLE    XI 
April    Lifting    Condensation   Levels 


E2IZ 

GMT 

OBS 

FC3T1 

FCST2 

FC5T3 

FC3TU 

13 

00 

705 

- 

- 

- 

- 

12 

3  70 

103 

- 

- 

- 

14 

00 

989 

59 

274 

- 

- 

12 

305 

- 

191 

919 

- 

15 

00 

506 

- 

- 

861 

234 

12 

a  t*  5 

— 

— 

— 

190 

E2TE 

GMT 

OBS 

FCST5 

FCST6 

FCST7 

FCST8 

15 

12 

445 

315 

- 

- 

- 

16 

00 

623 

423 

398 

- 

- 

12 

256 

- 

303 

297 

- 

17 

00 

- 

- 

- 

206 

78 

12 

- 

- 

- 

- 

0 

TABLE    XII 

April  Potential   Teaperatures 

DATS    GMT  OBS  FCST1  FCST2          FCST3          FCST4 

13  00  13.6  - 

12  7.3  12.8  - 

14  00  15.7  12.9  10.7 

12  8.2            -  12.5            15.1 

15  00  13.2           -  -  14.7            10.3 
12  10.6            -  -  -               11.4 

DATE    GMT  OBS  FCST5  FCST6          FCST7         FCST8 

15  12  10.6  9.9  - 

16  00  12.4  9.6  8.7 

12  8.3            -  7.1             12.3 

17  00                         -  -  9.8  5.6 
12  -  -  2.7 


Close  examination  reveals  that  the  potential  temperature 
forecasts  did  not  reflect  the  observed  diurnal  variation. 
The  12-hcur  forecast  generally  showed  cooling  or  heatir.q  as 
it  shculd  but  not  sufficiently  to  agree  with  cbserved 
results.  The  24-hour  forecast  tended  to  reverse  the  diurnal 
trend    and   to  have   coding  instead   of    heating    cr   vice-versa. 

It  is  believed  that  the  land  again  influenced  the 
cbserved  inversion  heights.  Figure  4.3  shows  the  cbserved 
and  fcrecasted  inversion  height  values  for  this  case.  Clcse 
examination   reveals      that   thare    appears      to    be      two    problems 
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TABLE    XIII 

April   Moisture   Yalues 


p  JfTT 

GMT 

OB5 

FCST  1 

1V 

00 

6.6 

- 

12 

5.2 

8.9 

m 

00 

6.5 

9.0 

12 

5.7 

- 

15 

00 

7.2 

- 

12 

6.  2 

— 

ESTE 

GMT 

OBS 

FCST5 

15 

12 

6.2 

6. a 

16 

00 

6.4 

5.9 

12 

5.9 

- 

17 

00 

- 

- 

12 

- 

- 

FCST2 


6.9 
8.2 


FCST6 

5.6 
5.3 


FCST  3 


6.5 
6.6 


FCST7 


7.6 
6.7 


FCST4 


6.8 

7.6 

FC  S  T  8 


5.4 
5.2 


which  are  a  combination  of  the  problems  discussed  fcr  the 
summer  and  winter  cases.  The  first  parr  of  the  period  was 
strongly  influenced  by  diurnal  fluctuations  and  the  fore- 
casted values  were  siiilar  to  these  in  the  summer  case.  The 
winter  case  was  influenced  by  subsidence  and  surface  winds. 
In  this  case,  it  is  believed  the  northerly  winds  caused  the 
increase  in  the  observed  inversion  at  15/0000  GMT  due  to 
tcpcgraphic  induced  large  mechanical  mixing.  The  vertical 
velocity  allowed  the  model  to  decrease  the  heights  fcr  that 
time. 

E.       DISCDSSION    OF    HODEI    RESULTS    FOR    PERIOD    17     (FALL) 

Tables  XIV  -  XVII  list  the  observed  and  forecasted 
values.  As  before  predicted  inversion  heights  for  this  case 
appear  tc  be  no  better  than  those  for  previous  cases.  The 
formation  of  clouds  was  forecasted  accurately  since  "the 
nodal  predicted  cloud  formation  at  a  time  when  low  clouds 
occurred  in  the  area.  However,  forecasted  lifting  condensa- 
tion   levels   were   not    in   agreement   with   those    observed. 
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TABLE    XI? 

August-September    Inversion    Results 

EATI    GMT    DBS       FCST1       RMS       FCST2       RMS       FCST3       RMS 


29 

CO 

287 

- 

- 

- 

- 

- 

- 

12 

297 

137 

16  0.0 

- 

- 

- 

- 

30 

00 

365 

67 

336.2 

524 

159.0 

- 

- 

12 

31  1 

- 

- 

638 

363.6 

3  33 

72.0 

31 

00 

391 

— 

— 

— 

— 

4  96 

127.3 

DATE 

GET 

OBS 

T?ccT4 

RMS 

FCST5 

RMS 

FCST6 

RMS 

31 

00 

391 

*  512 

121.0 

- 

- 

- 

- 

12 

553 

666 

165.6 

574 

21.0 

- 

- 

1 

00 

322 

- 

- 

755 

433.5 

819 

497.0 

12 

H 

- 

— 

— 

- 

1057 

- 

EATS 

GMT 

OBS 

FCST7 

RMS 

FCST8 

RMS 

FCST9 

RMS 

1 

12 

H 

363 

- 

— 

— 

- 

— 

2 

00 

356 

545 

18S.0 

- 

- 

- 

- 

12 

356 

- 

- 

273 

83.0 

- 

- 

3 

00 

4  1  1 

- 

- 

225 

244.2 

449 

38.0 

12 

469 

- 

— 

— 

5  09 

55.2 

EATE 

GMT 

OBS 

FCST1C 

)    RMS 

FCST11 

RMS 

FCST12 

RMS 

3 

12 

469 

491 

22.0 

— 

- 

- 

- 

a 

00 

422 

714 

292.8 

684 

262.0 

- 

- 

12 

529 

- 

— 

8  54 

417.5 

454 

75.0 

5 

00 

424 

— 

— 

— 

— 

562 

157.1 

DATE 

GMT 

OBS 

FCST13 

I    RMS 

FCST14 

•    RMS 

FCST15 

5 

00 

424 

■     675 

2510 

- 

- 

- 

12 

429 

823 

467.; 

2    404 

25 

- 

6 

00 

- 

- 

- 

496 

- 

4  98 

12 

t 

- 

- 

- 

- 

561 

The  forecasted  specific  humidity  values  have  an  inter- 
esting trend  (shewn  in  Table  XVII) .  Forecasts  initialized 
from  evening  data  *ere  good  but  those  initialized  from 
morning  data  were  bad.  When  initialized  with  evening  data, 
the  predicted  specific  humidity  values  were  representative 
cf      the   ctserved     values.  Obviously,         pcor   forecasts     of 

specific  huiridity  or  any  other  key  parameter  lead  to  pcor 
prediction   of  inversion    heights. 

Pctertial  temperatures,  in  all  but  two  runs,  were  fore- 
cast to  always  decrease  (cocl)  with  time.  The  two  morning 
to  evening  forecast  cases  did  not  reflect  sufficient  daytime 
heating,      even    with      clouds    present.         To    evaluate      the    pcor 


TABLE    XV 

august -Sap taster   Lifting   Condensation   Levels 

EATE    GHT  OBS  FCST1  FCST2         FCST3         FCST4 

29         00  731  - 

12  196  503  - 

2C         00  520  234  227 

12  206  -  0                 78                 - 

31         00  406  0  297 

12  364  -  206 

EATE  GMT  OBS  FCST5  ECST6    FCST7    FCST8 

2  1    12  36  4  105  ... 

1  00  467  227  613 

12  a  -  580      101 

2  00  331  -  150 

12  205  -  0 

3  00  577  -  0 

DATE  GMT  OBS  FCST9  ECST10   FCST11   FCST12 

3  00  577  237  - 
12  284  122  325 

4  00  496  -  452      424 

12  452  -  341        0 

5  00  525  -  0 

DATE  GMT  OBS  FCST13  ECST14   FCST15 

5  00  525  374 

12  310  263  394 

6  00  -  186  194 
12  -  0 


performance  in  potential  temperature  forecasts,  the  subsi- 
dence ratss,  specific  humidity  and  diurnal  radiaticn  varia- 
tions were  compared.  When  the  forecast  was  for  a  mcrning  to 
evening  period,  excluding  the  two  previously  mentioned  fore- 
casts, ccntinued  coding  was  predicted  when  heating  should 
have  been  predicted.  This  result  was  independent  of  subsi- 
dence rates  or  presence  of  clcuds.  Therefore,  for  mcrning 
to  evening  forecasts,  radiation  effects  were  not  sufficient 
to  introduce  the  needed  daytime  heating.  From  evenina  to 
morning      the  major      problem      was      too    much      cooling.  This 

happened  in  five  of  the  forecasts  with  another  five  fore- 
casts shewing  goGd  ccoling  rates  but  bad  forecasts  for  the 
24-hcur  point  since  the  previous  12-hour  forecast  was  bad. 
The  geed  ceding  rates  occurred  when  the  subsidence  was  at 
its  largest  compared  to  the  excessive  lowering  forecasts. 
This      would     indicate      possible    problems      in      the      radiation 
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TABLE    XVI 
August-September    Potential  Temperatures 


EATE 

GMT 

OBS 

FCST1 

FCST2 

FCST3 

FCSTU 

29 

00 

17.  8 

- 

- 

- 

- 

12 

13.  5 

16.6 

- 

- 

- 

30 

00 

16.  6 

15.5 

10.4 

- 

- 

12 

13.8 

- 

7.6 

13.7 

- 

21 

00 

15.  e 

- 

- 

11.5 

9.3 

12 

14.2 

— 

— 

— 

7.6 

EflTE 

GMT 

OBS 

FCST5 

FCST6 

FCST7 

FCST8 

21 

12 

14.  2 

12.3 

- 

- 

- 

1 

00 

17.  6 

13.2 

12.1 

- 

- 

12 

M 

- 

10.3 

13.8 

- 

00 

16.  1 

- 

- 

12.5 

— 

12 

12.8 

- 

- 

- 

8 .  6 

3 

00 

15.  4 

— 

— 

— 

9.7 

DATE 

GMT 

OBS 

FCST9 

FCST10 

FCST11 

FCST12 

3 

00 

15. a 

9.9 

- 

- 

- 

12 

13.  1 

7.8 

12.2 

- 

— 

4 

00 

16.  C 

- 

10.8 

10.1 

- 

12 

14.  1 

- 

- 

7.9 

12.2 

5 

00 

16.  6 

— 

- 

— 

11.2 

DATE 

GMT 

OBS 

FCST1 3 

FCST14 

FCST15 

5 

00 

16.  € 

11.6 

- 

- 

12 

13.8 

9.7 

16.1 

- 

6 

00 

- 

- 

12.8 

11.2 

12 

- 

- 

- 

8.7 

segment  which  wculd  be  cooling  the  potential  temperatures 
for   this   tiire   of   day. 

Figure  4.4  shows  the  observed  and  forecasted  inversicns 
for  this  period.  Fey  factors  affecting  inversion  heights 
during  this  period  were  diurnal  heating  variations,  surface 
winds,        and  clouds.  The    temperatures      showed    a      definite 

diurnal  pattern  throughout  the  period  as  did  the  lifting 
condensation  level.  The  inversion  displayed  a  diurnal  trend 
only  part  of  the  time,  29/1200  GMT  through  31/0000  GMT  and 
3/0000   GMT    through    the   remainder   of   the    period. 

The  low  subsidence  values  contributed  greatly  to  the 
high  inversion  height  values  forecasted  during  this  period. 
Of  the  remaining  forecasts,  those  with  heights  forecast  to 
decrease  were  evening  to  morning  forecasts  with  subsidence 
occurring.  The  observed  heights  increased  (except  for  the 
2/O00C   GKT      value    which      remained   constant)         due    to      strong 
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TABLE  XVII 

August-September  Moisture  Values 


DATE  GMT 
2S    00 

12 
3C    00 

12 
31    00 

12 

DATE  GMT 
31    12 

1  00 
12 

2  00 
12 

3  00 

EiTE  GMT 
3    00 

12 
I*    00 

12 
5    00 

DATE  GMT 

5  00 
12 

6  00 
12 


OBS 
8.5 
3.7 
3.  1 
8.8 
9.0 
8.  3 

OBS 
8.3 
9.9 

M 
9.6 
8.2 
8.0 

OBS 
8.0 
8.  1 
8.7 
7.9 
9.0 

OBS 
9.0 
8.  3 


FCST1 

9.0 
9.7 


FCST5 

8.  4 
8.3 


FCST9 
6.7 

6.1 


FCST1 3 
6.9 
6.4 


FCST2 


6.9 
6.6 


FCST6 

6.3 
5.7 


7.4 
6.3 


9.3 

8.3 


FCST3 


9.3 
9.3 


FCST7 


9.3 
8.3 


6.  1 
5.5 


FCST14   FCST15 


FCST4 


6.4 
5.8 

FCST8 


9.  1 
8.9 


FCST10   FCST11   FCST12 


9.0 
9.  1 


7.5 
7.0 


winds  (weaker  and  more  westerly  at  2/0000  GMT)  which  lifted 
the    inversion   due   to    nixing. 

At  31/1200  GMT  an  observed  height  decrease  occurred  for 
the  morning  to  evering  forecast  time  but  the  forecasts 
showed  rises.  Tc  explain  this  requires  consideration  of  all 
three    factors   mentioned   above.  The    observed   height    values 

decreased  dua  to  the  increase  in  cloud  coverage  and 
continued  onshore  flew  whereas  the  normal  situation  is 
effshcre  flew  due  to  diurnal  changes.  The  radiation  compu- 
tations also  appear  to  have  been  responsible  for  errors  in 
the    forecasts   for   this   period.  Radiation   was   not    the   only 

problem  since  erroneous  specific  humidity  and  vertical 
velocities   contributed.  There    was      land   influence      on  the 

parameters  involved  in  height  predictions,  as  in  the  ether 
cases  . 
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Figure  4. 4    Graphical  Depiction  of  Fall  Case 
Observed  ana  Forecasted  Inversion  Heights. 
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I.        SOMMJBI    AND    GENEEfil    COSBENTS 

It    is    easily      seen   by   the    large    root      mean   square    secies 

a  r\  A  *  Y  o     Is  res  ^i'  f  ■*  csprog         fo  S  t  W  S  S  n  forecast     3.  Ti  d        CTS'.tr  V°d 

trends,  depicted  graphically,  that  the  model  did  net  yield  a 
good  pradiction  cf  inversion  heights  at  Vandenberg  AFB.  The 
model  also  did  not  yield  very  good  forecasts  cf  the  ether 
paraceteis. 

The  discussions  present  an  explanation  of  the  prckatle 
components  where  errcrs  occurred.  There  are  certain  prime 
factcis  v»hich  cause  the  greatest  problem.  The  largest  error 
producer  is  the  fact  that  Vanbenberg  AFB  is  a  land  base 
which  strcngly  experiences  diurnal  variations  in  temperature 
and  winds.  The  primary  predictive  values  of  the  model 
(subsidence,  entrain  snent,  radiation  and  presence  of  clcuds) 
are  also  influenced  ty  the  land  effect.  Daytime  heating 
which  increases  temperature  and,  hence,  dissipates  clouds  is 
not  accounted  for  in  the  MAEL  model. 

Topography  and  certain  wind  directiens  increase  the 
turbulence  in  the  mixed  layer.  With  mountains  north  as  well 
as  south  of  Vandenberg  AFB,  this  condition  can  result  when 
the  effshcre  flow  dees  not  set  up  during  the  night  and 
strong  northerly  or  southerly  flew  continues.  This  leads  to 
deepening  cf  the  mixed  layer  and  can  overcome  the  effects  of 
subsidence.  The  result  of  this  would  be  to  lift  instead  of 
lower  the  inversion. 

It  is  concluded  that  this  study  yielded  conclusive 
evidence  that  the  MAEI  model  is  not  appropriate  for  use  at  a 
coastal  station  because,  even  though  coastal,  land  influ- 
ences play  a  significant  role.  The  results  from  the  study 
ty  Gleascn  (1982)  reveals  definite  contrasts  with  this  one. 
He  had  geed  results  for  San  Nicolas  Island  and  was  able  to 
show  that  the  Q-methcd  of  calculating  vertical  velocities 
was  test. 
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v-     SOHSARX    ASD    CONCLUSIONS 

Accurate  predict ions  of  variations  in  the  height  of 
atmospheric  inversions  are  a  necessity  in  proper  usage  of 
modern  weaponry.  This  is  because  of  the  impact  changes  in 
inversion  heights  have  on  the  ducting  of  weapons  guidance 
system  or  of  tracking  system  signals.  Atmospheric  parame- 
ters are  continually  changing  and  those  parameters  which 
govern  inversions  can  easily  change  their  entire  profile  in 
a  matter  cf  hours.  forecasts  of  changes  in  inversion  prop- 
erties irust  start  «ith  forecasts  of  those  parameters, 
namely,  temperature  and  specific  humidity.  To  do  this  the 
narire  atmospheric  boundary  layer  model  used  at  the  Naval 
Postgraduate  School    was   utilized. 

Data  periods  which  covered  each  season  were  chosen  for 
the  Vandentsrg  AF3  area.  All  four  cases  were  dominated  by  a 
high  pressure  system  over  the  eastern  North  Pacific  Ocean. 
Clear  as  well  as  clcudy  days  were  considered  to  allow  full 
range  testing  for  a  station  strongly  influenced  by  land 
conditions.  Data  for  these  periods  were  from  rawinsonde 
observations,  satellite  information,  and  surface  and  500  mb 
synoptic   analyses.  Model    input      parameters   were      from   the 

rawinscr.de  observations  and  satellite  derived  sea-surface 
temperatures. 

Results  from  each  model  run  were  compared  with  observed 
values  to  determine  accuracies  of  the  model.  It  was  found 
that  the  model  did  net  yield  accurate  predictions  during 
most  cf  the  cases.  Final  interpretations  of  these  inaccura- 
cies always  point  to  strong  land  influence.  Diurnal  varia- 
tions influenced  by  the  land  are  much  greater  than 
variations  over  water  and  are  not  reflected  in  the  model 
calculations. 
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The  end  result  cf  this  study  is  that  this  MAEL  modal 
dees  net  work  for  Vandenberg  afb.  It  is  not  currently  valid 
for  staticrs  strongly  influenced  by  land  induced  diurnal 
variations.  The  redevelopment  cf  the  model  to  treat  such 
factors  is  required  tc  properly  handle  the  land  influences. 
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